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Abstract. Observations of the shell-type supernova remnant SN1006 have been carried out with the H.E.S.S. 
system of Cherenkov telescopes during 2003 (18. 2h with two operating telescopes) and 2004 (6.3h with all four 
telescopes) . No evidence for TeV 7-ray emission from any compact or extended region associated with the remnant 
is seen and resulting upper limits at the 99.9% confidence level are up to a factor 10 lower than previously- 
published fluxes from CANGAROO. For SN1006 at its current epoch of evolution we give limits for a number 
of important global parameters. Upper limits on the 7-ray luminosity (for i5=0.26 to 10 TeV, distance d=2 kpc) 
of L 7 < 1.7 x 10 33 erg s _1 , and the total energy in corresponding accelerated protons, W p < 1.6 x 10 50 erg (for 
proton energies E p ~ 1.5 to 60 TeV and assuming the lowest value n = 0.05 cm -3 of the ambient target density 
discussed in literature) are estimated. Extending this estimate to cover the range of proton energies observed in 
the cosmic ray spectrum up to the knee (we take here E v ~ 1 GeV to 3 PeV, assuming a differential particle 
index —2) gives W p < 6.3 x 10 50 erg. A lower limit on the post-shock magnetic field of B > 25/iG results when 
considering the synchrotron/inverse-Compton framework for the observed X-ray flux and 7-ray upper limits. 
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1. Introduction 

Observations of shell-type supernova remnants (SNR) at 
multi-GeV to TeV 7-ray energies have long been mo- 
tivated by the idea that they are the prime source of 
hadronic cosmic-ray (CR) acceleration in our galaxy (see 
for example Drury, Aharonian & Volk 119941 and Naito & 
Takahara 1994 ). The detection of non-thermal X-ray emis- 
sion from the young shell-type SNR SN1006 (Koyama et 
al. Allen et al. r27)0l")> suggests that SN1006 is a site 

of electron acceleration to multi-TeV energies. Subsequent 
arc-second resolution results from Chandra (Long et al. 
2003 , Bamba et al. 12003} reveal the presence of several 
bright non-thermal X-ray arcs concentrating in the NE 
and SW regions of the SNR, which likely trace the location 
of strong shocks where at least electrons are accelerated. 
These, and similar X-ray results from a number of other 
SNR over the past ten years, have provided a priority list 
of such targets for ground-based 7-ray telescopes. The re- 
sults of CANGAROO-I (Tanimori et al.QSSHJ, suggesting 
TeV 7-ray emission from the NE rim of SN1006 supported 
the notion that this SNR is capable of electron accelera- 
tion up to energies ~ 100 TeV, when interpreted in the 
synchrotron/inverse-Compton (IC) framework. The TeV 
photon spectrum is compatible with a power law of pho- 
ton index T -2.3 (for dN/dE ~ £~ r ) in the energy range 
1.5 to 20 TeV (Naito et al. IT0001 Tanimori et al. 1277711)) . 
Later data from the CANGAROO-II 10 metre telescope 
also revealed SN1006 as an emitter of TeV 7-rays (Hara 
et al. 12001(1 . yielding a compatible energy spectrum. The 
stand-alone HEGRA CT1 telescope has also reported a 
significant excess in observations taken at very large zenith 
angles and therefore at a high energy threshold E > 18 
TeV (Vitale et al. 12005)1 . 

The interpretation of these X-ray and ground-based 
7-ray results generally involves the electronic syn- 
chrotron/IC and/or hadronic 7r°-decay channels. Either 
channel can dominate, depending on, for example, the 
density of ambient matter n and magnetic field strength 
B. Detailed sampling of the 7-ray spectra and morphol- 
ogy are required to disentangle the electronic and hadronic 
components. The H.E.S.S. (High Energy Stereoscopic 
System) experiment, the first of the next-generation 
ground-based 7-ray detectors to utilise the stereoscopic 
technique, has observed SN1006 during 2003 and 2004. 
The prime motivation of these observations, apart from 
confirmation of the CANGAROO results, has been to de- 
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Table 1. Summary of H.E.S.S. observations of SN1006. 



Obs. 


Runs 


ON source 


a < z > 


Telescopes 


Period 




livetime [h] 




in use 


Mar 2003 


18 


4.1 


23.9 


2 


Apr 2003 


21 


5.0 


25.4 


2 


May 2003 


25 


9.1 


23.3 


2 


2003 Total 


65 


18.2 


24.0 




May 2004 


15 


6.3 


28.5 


4 


All Total 


80 


24.5 


25.0 




a: Average zenith 


angle [deg 


;] of observations. 







termine the nature of any 100 GeV to TeV 7-ray emission 
from SN1006. 

Operating in the Southern Hemisphere, the H.E.S.S. 
experiment consists of four identical Cherenkov telescopes 
each with mirror area ~107 m 2 (Bcrnlohr et al. [2003 
Cornils et al. 120051 Vincent et al. 12005)) . All four tele- 
scopes have been operating since December 2003. Here, 
we present results using H.E.S.S. data taken during 2003 
(two telescopes operating in stereo but with considerable 
dead time of ~50% due to the lack of a central trigger sys- 
tem EES). The four-telescope mode of 2004 with the cen- 
tral trigger is the most sensitive configuration for H.E.S.S. 
We update here previous results of Masterson et al. ( 2003) 
which utilised a part of the 2003 dataset analysed without 
stereoscopic reconstruction (i.e. mono-mode). The large 
field of view (FoV ~ 5° diameter) of H.E.S.S. easily en- 
compasses the SN1006 shell, and the angular resolution 
of ~ 0.1° attained (event-by-event) permits a search for 
7-ray sources from different regions associated with the 
SNR. 

2. Analysis & Results 

H.E.S.S. observation runs of ^28 min duration were taken 
using the so-called wobble mode in which the tracking po- 
sition is displaced ±0.5° in declination with respect to 
the centre of SN1006 (RA 15 /l 02 m 48.4 s Dec -41°54'42" 
J2000.0). Runs were accepted for analysis if they met a 
number of quality control criteria based on the recorded 
CR rate, the number of malfunctioning pixels in each cam- 
era and also the tracking performance. In particular tele- 
scope tracking was accurate to ~20 arcsec for these data 
and verified by comparing the locations of bright stars in 
the camera fields of view. Table ^ summarises observa- 
tions for those runs meeting quality criteria. Without cor- 
rection for system deadtime the overall observation times 
were 30h (2003, with - 50% deadtime) and 7h (2004, with 
- 10% deadtime). A total of 18. 2h (2003) and 6.3h (2004) 
livetime of ON source data (where SN1006 is in the FoV) 
were available for analysis. 

Individual telescope events coincident in time were 
merged for stereo analysis, followed by Cherenkov im- 
age reduction. Image reduction employs image 'cleaning' 
which removes camera pixels dominated by skynoise, flat 
fielding of the camera responses using a LED flasher, and 
absolute conversion from pixel ADC counts to photoelec- 
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tron units using conversion constants obtained from spe- 
cial low-illumination runs. Image moments such as width 
and length according to Hillas 1|1985[1 are used as a ba- 
sis to reject the dominating CR background. Rejection 
of the CR background is achieved by application of im- 
age shape cuts [mean-reduced- scaled-width MRSW and 
mean-reduced-scaled-length MRSL). A further cut, 0, the 
difference in assumed and reconstructed event directions is 
also applied. Event directions are reconstructed according 
to algorithm '1' as detailed by Hofmann et al. I|1999jl . For 
a point-like source the combination of cuts on shape, and 
direction cut < 0.14° allows for rejection of over 99.9% 
of CRs whilst retaining ^40% of 7-rays above a threshold 
energy of 260 GeV (2003) and 110 GeV (2004). The full 
data processing chain and CR rejection applied to these 
data were defined a-priori using Monte Carlo simulations 
of 7-ray and real CR events, and fully verified on the Crab 
Nebula. Details can be found in Aharonian et al. ((2004). 

The skymap in Fig. ^ presents the excess significance 
for 7-ray like events (i.e. after cuts) over the RA/Dec plane 
centred on SN1006. Data from 2003 and 2004 were com- 
bined for these analyses. At each bin position, 7-ray-like 
events are summed within a circle of radius cu t < 0.14° 
and the CR background is estimated from a ring region 
of radius 0.5° surrounding the source region. The ring 
region is chosen to give a solid angle ratio a of 1:7 be- 
tween the source and background regions respectively. 
The skymap bins of Fig. ^ ar e therefore correlated. The 
cut 9 cu t < 0.14° is appropriate for point-like 7-ray and 
marginally extended (~ few arcmin) emission according to 
the H.E.S.S. point spread function (PSF). No significant 
excess suggesting 7-ray emission is evident in the skymap. 
A very similar skymap is obtained when using the alter- 
native template model ; Howell 2003) as a CR background 
estimate. Event statistics and 7-ray flux upper limits at 
the 99.9% confidence level, summarised in Tabled have 
been calculated for a number of a-priori chosen locations 
based on the X-ray morphology as imaged by ASCA and 
the results of CANGAROO and HEGRA CT1. The up- 
per limits in the table have assumed a differential photon 
index of T=2, and later (in Fig. EJ) we give the upper limit 
band for a range of assumed indices (2 to 3) . The a-priori 
locations, also indicated in Fig. ^ are defined as: 

1. H.E.S.S. Point: Upper limit on point-like emission 
(E > 0.26 TeV) at the CANGAROO position using 
9 cut < 0.14° (appropriate for the H.E.S.S. PSF). This 
is also used as an upper limit for the NE Rim position. 

2. CANGAROO Point: Upper limit at the 
CANGAROO position using a cut < 0.25° cut 
appropriate for the CANGAROO 3.8m PSF and 
selecting events above their energy threshold [E > 1.7 
TeV). 

3. CANGAROO Point 2003: As above using only 

2003 data. 

4. CANGAROO Point 2004: As above using only 

2004 data. 



5. High Energy Point: Excess events at the 
CANGAROO position using cut < 0.14° cut appro- 
priate to the H.E.S.S. PSF but with an energy thresh- 
old E > 18 TeV to compare with the HEGRA CT1 re- 
sult (Vitale et al. l2003|l . Due to low statistics, the back- 
ground is estimated from the full FoV (radius 2.0°) ex- 
cluding the source region and corrected for differences 
in expected FoV efficiency. The dominant error comes 
from the on count statistics. 

6. Whole SNR: Entire SNR at the remnant centre us- 
ing cut < 0.50°, and H.E.S.S. energy threshold of 
(E >0.26 TeV). 

7. SW Rim: As for position 1 but centered on the SE 
rim. 

3. Discussion & Conclusions 

Comparing the H.E.S.S. upper limits with the fluxes from 
CANGAROO-I (1996 and 1997 observations) reveals a 
discrepancy of a factor ~10 (see Fig. EJ). It is also clear 
that no 7-ray emission arises in just one of the years that 
H.E.S.S. has observed, 2003 nor 2004 (see results given 
in Table 12 for the CANGAROO position, 2003 and 2004 
respectively). Thus H.E.S.S. observations do not confirm 
the previously-published CANGAROO fluxes for SN1006. 
One would have to invoke a 7-ray flux variation over 
timescales less than 1 percent of the SNR age of an or- 
der of magnitude, in order to explain the non-detection 
by H.E.S.S. This is however quite unlikely to occur within 
the long-accepted diffusive shock acceleration framework 
for particle acceleration in SNR. 

The HEGRA CT1 telescope observed SN1006 at very 
large zenith angles and obtained a preliminary flux of F 
(E> 18 TeV) = 2.5±0.5 sta txl0" 13 ph cm^s" 1 (Vitale et 
al. I2003[) . We applied a matching high threshold cut to 
H.E.S.S. data and obtained zero counts on source, (with 
s — ab=— 0.14 excess events). The 99.9% upper limit on 
this is a Poisson distribution with a mean of 7 events 
(we take here the dominant error in the on source counts 
s and quote the Poisson error). In comparing with the 
HEGRA CT1 result one must also take into considera- 
tion the systematic uncertainty on this flux, estimated to 
be >35% (Vitale et al. 2003). Given our exposure and 
considering both statistical and (uncorrelated) systematic 
uncertainties, the HEGRA CT1 flux would yield 11±6 ex- 
cess counts if seen by H.E.S.S. The lower fluctuated value, 
at 5 counts, is consistent with our observed zero counts 
with a chance probability of ~ 7 x 10~ 3 . We also include 
in Fig. H the 99.9% upper limit for the whole SNR as a 
function of energy. This allows comparison in cases where 
emission from models is predicted from regions not neces- 
sarily restricted to the NE and SW rims. 

The H.E.S.S. upper limits permit estimation of im- 
portant parameters for SN1006, in the framework of this 
source as a particle accelerator. One can firstly estimate 
an upper limit to the 7-ray luminosity as L~ ( = AirFd 2 
erg s _1 for a 7-ray flux F and source distance d. Distance 
estimates for SN1006 are in the range d ~0.7 to 2.0 kpc, 
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Fig. 1. Skymap of excess significance (according to Eq. 17 of Li & Ma (1983)) centred on SN1006. Indicated are 
all a-priori integration regions used in these analyses (numbered according to the text and Table EJ) and ASCA X- 
ray contours highlighting the X-ray emission from the NE and SW rims. The distribution of significances is fit by a 
Gaussian distribution with a mean n — —0.016 ± 0.004 and sigma a = 1.035 ± 0.003. 



(see Allen et al. 2001 for a recent summary of measure- 
ments). Taking the H.E.S.S. upper limit for the whole SNR 
(E > 0.26 TeV), distance d = 2.0 kpc, assuming a differ- 
ential photon spectrum of T — 2 and energy range 0.26 
to 10 TeV gives L 7 < 1.7 x 10 33 erg s -1 . This applies to 
the present epoch of SNR evolution for SN1006 which is 
likely in the early Sedov phase (Berezhko et al. 12002(1 . We 
also estimate the energy in the corresponding accelerated 
protons W p over the SNR lifetime as W p ~ L 7 r pp erg. The 
W p estimate applies to protons of energy E p ~ 1.5 to 60 
TeV, and accounts for the average inelasticity (~ 0.17) of 
converting energy from protons to 7-rays in this scheme. 
The rather energy independent characteristic lifetime r pp 
of accelerated protons interacting with an ambient mat- 
ter of density n cm~ 3 (p + p — > tt° + X — > 2j + X) is 
Tp P ~ 4.5 x 10 15 (n/cm -3 ) -1 sec. With L 1 derived above 
we arrive at W p < 7. 8x10 (n/cm -3 ) -1 erg. Values for 
n are in the range n — 0.05 to ~ 0.3 cm -3 based on X- 
ray and optical observations. Adopting the lowest value 
n = 0.05 cm' 3 of the range of densities yields the upper 



limit W p < 1.6 x 10 50 erg. 



The choice of n in the W p esti- 



mate should reflect the average density of ambient matter 



inside the remnant undergoing interaction with the down- 
stream cosmic-rays (which are assumed to be confined in- 
side the SNR). A lower limit on this value would be that 
of the ambient matter. We can account for the fact that 
the energy band of H.E.S.S. observations corresponding to 
protons (~ 1.5 to 60 TeV) is a fraction ~ 0.25 of the total 
expected energy in protons over a much wider range of 
energies covering the cosmic-ray spectrum up to the knee 
(we assume here ~ 1 GeV to 3 PeV with a differential 
particle index of —2). Given this, an upper limit on the 
full proton budget would be W p < 6.3 x 10 50 erg. 

Under the synchrotron/IC (scattering on the ubiqui- 
tous cosmic microwave background) scenario, knowledge 
of the synchrotron X-ray flux f x and TeV IC flux / 7 aris- 
ing from the same electrons permits a clear constraint on 
the B field in the shock-compressed regions (downstream 
from the shock) of the SNR according to f 1 {E 1 )/f x {E x ) ~ 
0.1(B/10^G)~ 2 £ (Aharonian et al. ITW7|l . It is assumed 
here that the X-ray and 7-ray emitting regions are of 
the same size, and thus the filling factor £ ~ 1. Such 
a case arises for high B fields when electrons rapidly 
cool at their synchrotron production regions. High val- 
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Table 2. Event statistics and flux upper limits (using the method of Feldman & Cousins 1998, assuming a differential 
photon index F" = 2) for a number of a-priori-chosen regions associated with SN 1006: s - source events; b - CR 
background estimate; Normalisation factor a; See text for an explanation of the various regions. Except for limits 3 
and 4, numbers are derived exclusively from combined 2003 and 2004 H.E.S.S. data. 



Region 


a RA 


a Dec 


Pcut 


s 


b 


a 


b S 


c ,99.9% 


1. H.E.S.S. Point (E > 0.26 TeV) 


15 h 03 m 48" 


-41°45' 


0.14° 


1072 


6854 


0.150 


1.19 


2.19 


(also NE rim) 


















2. CANGAROO Point (E > 1.7 TeV) 


15 h 03 m 48 3 


-41°45' 


0.25° 


127 


792 


0.153 


0.48 


0.34 


3. CANG. Point (E > 1.7 TeV) 


15 h 03 m 48 s 


-41°45' 


0.25° 


56 


309 


0.152 


1.2 


0.45 


(2003 H.E.S.S. Data) 


















4. CANG. Point (E > 1.7 TeV) 


15 h 03 m 48 s 


-41°45' 


0.25° 


71 


483 


0.154 


-0.35 


0.49 


(2004 H.E.S.S. Data) 


















5. H.E. Point (E > 18 TeV) 


15 h 03 m 48 3 


-41°45' 


0.25° 





20 


0.0069 




0.18 d 


6. Whole SNR (E > 0.26 TeV) 


15 h 02 m 48.4 s 


-41°54'42" 


0.50° 


13358 


63113 


0.217 


-2.8 


2.39 


7. SW Rim (E > 0.26 TeV) 


15 h 02 m 4 s 


-42°06'3" 


0.14° 


1042 


6754 


0.150 


0.80 


1.98 



a. RA & Dec J2000.0 epoch 

b. Significance from Eq. 17 of Li &c Ma^^ using normalisation factor a 

c. 0^' 9% = 99.9% integral upper limit (xlO -12 ph cm _2 s _1 ) 



d. Dominant error is from s. Wc quote here the 99.9% upper limit, which for a Poisson of mean zero is 7 counts. 



ues for B > 40/iG are in fact implied by X-ray ob- 
servations and subsequent interpretation in the diffusive 
shock acceleration framework (Allen et al. 120011 Berezhko 
et al. l2"no"2l IM)3l Bamba et al. IM)3l Yamazaki et al. 
I2UU1 Ksenofontov et al. I2UU3|) . With the H.E.S.S. up- 
per limits on the 7-ray emission, one can therefore es- 
timate a lower limit on the _B-field. The synchrotron 
(E x ) and IC (E^) photon energies are coupled accord- 
ing to E 1 ~ 1. 5{E x /0. lkeV)(S/10^G)" 1 TeV, requiring 
that the fluxes f x and / 7 in appropriate energy ranges 
be compared. Comparing the available X-ray energy flux 
(E = 0.1 to 2 keV) f x = 1.42 x 10~ 10 erg cm^s- 1 
(Allen et al-EDOU with that of H.E.S.S. ( upper limit en- 
ergy flux for the whole SNR over an appropriate energy 
range E~l to 10 TeV and assuming a —2.0 spectral in- 
dex) at / 7 =2.29x 10 -12 erg cm _2 s~ 1 yields a lower limit of 
B > 25/iG. This lower limit on B is consistent with values 
discussed earlier which result from comparisons of X-ray 
observations with detailed theory. In summary, H.E.S.S. 
observations of SN1006 have not revealed evidence for TeV 
7-ray emission. The resulting upper limits will be valuable 
in constraining further the parameters of 7-ray production 
in SN1006 (see e.g. Aharonian & Atovan 119991 Berezhko 
et al.|MEl Ksenofontov et al. lWRjl . 
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Fig. 2. Top: Comparison of experimental results for 
SN1006 from H.E.S.S. (integral upper limit 99.9% c.l. 
band at the CANGAROO position using 6> cut < 0.23°, see 
Table 0), CANGAROO-I fluxes of Tanimori et al. ijT^I 
and the HEGRA CT1 flux of Vitale et al. (IStMll. The solid 



line joins the two CANGAROO-I fluxes and is extended 
up to the HEGRA CT1 flux. The high energy H.E.S.S. up- 
per limit is shown for comparison with the HEGRA CT1 
flux. The H.E.S.S. UL band arises when assuming a dif- 
ferential spectral index from T=2 to 3. Bottom: H.E.S.S. 
upper limit band for the whole SNR, after cuts described 
in the text. The band arises when assuming a differential 
spectral index of T—2 to 3. 
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